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Abstract: We study the observability of the Higgs boson in the "charming Higgs" model. 

In this model the Higgs boson primarily undergoes a cascade decay to four charm quarks via 
light intermediate pseudoscalars. Such a decay allows the Higgs boson to escape the most 
stringent LEP bounds on the Standard Model Higgs boson mass. If the light pseudoscalars 
are sufficiently light they become highly boosted and their decay products collimated into 
jets. We show that by using jet substructure techniques, this model is potentially observable 
at the LHC. For a Higgs boson mass of 100 GeV and light pseudoscalar mass of 12 GeV, we 
find a signal significance of 3.8<r with a luminosity of 30 fb _1 and that a 5<r significance can 
r\ . be obtained with 50 fb _1 of luminosity at the 14 TeV LHC. 
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1. Introduction 

The mechanism of electroweak symmetry breaking (EWSB) is one of the central unanswered 
questions of particle physics. While the Standard Model (SM) with one Higgs doublet provides 
a simple EWSB mechanism that is so far in agreement with experiment, LEP has already 
excluded a SM Higgs with mass below 114.4 GeV 0. Also, the ATLAS and CMS experiments 
at the Large Hadron Collider (LHC) are quickly excluding the SM Higgs. At 95% confidence 
level, the ATLAS experiment currently rules out the mass ranges 110 — 117.5, 118.5 — 122.5, 
and 129 - 539 GeV g, and CMS excludes 129 - 525 GeV @. As the limits on the SM Higgs 
boson tighten, we are motivated to investigate other possible mechanisms for EWSB. Two 
such examples are supersymmetry and technicolor 0]. 

One may also note that electroweak precision constraints prefer a SM Higgs with mass at 
or slightly below 100 GeV ||, beneath the LEP direct detection bounds. Hence, a particularly 
interesting class of models are those that contain a SM-like Higgs boson that escapes the 
strictest LEP bounds and has mass nearer that preferred by precision measurements. This 
can be accomplished in theories with a SM-like Higgs boson that preferentially undergoes a 
cascade decay into light SM particles. In particular, if the Higgs boson dominantly decays 
via a cascade h — )■ 2rj — > 4X, where rj is a pseudoscalar and X are light SM particles, the 
LEP bounds can be significantly relaxed || [7J. 

There has been much recent work on detecting a Higgs boson decaying to two pseu- 
doscalars and the pseudoscalar subsequently decaying to 27, 2r, 2b, 2// and 2g |8|, ^, [l(J. If 
the pseudoscalars preferentially decay to two Vs, the Higgs mass bounds from LEP are only 



mildly reduced to 110 GeV [11|. Hence, to evade the most stringent LEP bounds, many 



models of this type require the mass of the pseudoscalar to be less than the bottom quark 



threshold, 2m&, eliminating the r\ — > bb decay. However, one example, the so-called "charming 
Higgs" model |1^] , has matter representations that suppress the couplings of the pseudoscalar 
to down type quarks and charged leptons. Consequently, the pseudoscalar decay to charm 
quarks is dominant over the decay to bottom quarks or r leptons, even for pseudoscalar 
masses above the bottom quark threshold. In such a scenario, the Higgs boson can still have 
a mass around mz and avoid the most stringent LEP bounds. 

In this paper we study the Higgs production in association with al^at the LHC with 
a subsequent Higgs decay h — > It] — >■ 4c. Specifically, we consider the case where the Higgs 
boson mass, m^, is below the LEP bounds and the 77 mass, m™, is above the bottom quark 
threshold, m v > 2m],. Hadronically decaying Higgs bosons are difficult to observe at the LHC 
due to the large QCD backgrounds. We therefore consider the scenario where m v <C rrih- In 
this case the 77s become highly boosted and their decay products are collimated into single 
jets. Our signal jets then originate from color singlet particles decaying to two quarks while 
the QCD background jets originate from massless colored partons. Hence, the signal jets have 
a mass scale and substructure distinct from the background. Using the techniques developed 
in Ref. |l3|], we decompose the jets into subjets and perform a substructure analysis. 

The paper is structured as follows. In Section ^ we outline the simplified model we utilize. 
We study the observability of the charming Higgs associated production with a vector boson 



via a complete signal and background analysis in Section [3|. In Section |3.l| , we demonstrate 
that a subjet analysis improves the observability of the charming Higgs at the LHC. Finally, 
in Section [| we summarize our results and conclude. 

2. Simplified Model 



In the models of interest |12|, 14|, the Higgs boson arises as a pseudo-Goldstone boson (pGB) of 
an approximate global symmetry. The symmetry breaking pattern of these models gaurantees 
an additional light SM singlet pseudoscalar, 77. The pseudoscalar then has derivative couplings 
to the SM-like Higgs boson, h: 

2 vew / -. v EW 



where vew = 174 GeV is the electroweak breaking scale, and / is a global symmetry breaking 
scale. For / values on the order of the electroweak breaking scale, h will dominantly decay 
to two 77's. For example, if / = 350 — 400 GeV and the Higgs boson mass, nth, is around the 
Z-mass, the branching ratio of h to 1r\ is 80 — 90% and h to bb is 10 — 20%, consistent with 



the LEP bounds [12, 14]. 



The couplings of the pseudoscalar to a SM fermion are of the form 

iVfVhsf- (2.2) 

Typically, 77 couples most strongly to top and bottom quarks. Hence, for m^ > 2m;,, the 
77 predominantly decays into bb and the LEP bound is only slightly alleviated to nih > 



110 GeV [plj| . However, in the so-called "charming Higgs" model [[l4|], the 77 coupling to 
bottom quarks is suppressed by higher order operators. Also, r/ couples to rs through the 
mixing of r with a heavy partner, whereas the coupling to charm quark does not suffer from 
these suppressions. Hence, the dominant decay mode is then 77 — > cc for all values of m^. 

For simplicity and clarity, throughout the rest of this paper we take the branching ratios 
of the Higgs boson BR(/i — > rjrj) = 1 and pseudoscalar BR(t7 — > cc) = 1. For smaller branching 
ratios, the signal cross sections can be simply scaled. 

3. Signal and Background 

In the context of the simplified model in Sec. 0, we study the production of a Higgs boson 
in association with a W boson, pp — > Wh, with the subsequent decay h — >■ 2r/ — )■ 4c. Addi- 
tionally, we consider events where the vector boson decays to leptons, which allows our signal 
events to be more easily triggered and avoids pure QCD backgrounds. We allow the W to 
decay to both electrons and muons. The associated production with a Z-boson is also possi- 
ble; however, the branching ratio of Z to leptons is less than W to leptons and so we focus on 
purely Wh production. Similar studies have been performed for h — > Ag and h — > 46 ||, |l0|| . 
All signal and background events are generated using MadGraph5 [15| and then showered 



using Herwig 6.5.10 [16|. The charming Higgs model is incorporated into MadGraph using 



FeynRules vl.6 [17]. Jets are clustered using Fast Jet v2.4.2 [18] and the kx algorithm [IS] 
with a radius R = 0.5. 

For the purposes of this study, we consider a Higgs mass of 100 GeV and an 77 mass 
of 12 GeV. With this parameter choice the 77s will be highly boosted and their individual 
decay products will be indistinguishable. Hence, our signal consists of a lepton, two jets, and 
missing energy. After showering and clustering, we require at least two jets to pass the basic 
acceptance cuts 

p> T > 30 GeV \y j \ < 2.5, (3.1) 

where pt and y are transverse momentum and rapidity, respectively. Additionally, to trigger 
on the signal we apply the lepton and missing transverse energy cuts 

p e T > 30 GeV l/l < 2.5 fr > 25 GeV. (3.2) 

The signal cross section after these cuts is shown in the second column of Table 0. 

In Fig. H we present the signal invariant mass distributions of (a) the dijet system consist- 
ing of the two hardest jets, rrijj, and (b) the individual jet masses, rrij, for the hardest (solid) 
and second hardest (dashed) jets. The cuts in Eqs. Q3.1| ) and ( |3.2| ) have been applied. As 
expected, the invariant mass of the dijets is peaked at the higgs mass m,h = 100 GeV. Also, 
since each jet originates from a massive particle, the individual jet masses peak at rrij = m v . 
The differences in the hardest and second hardest jet mass distributions can be understood 



Table 1: Cross section for signal (Wh) and background, signal to background ratio, and signal 
significance with consecutive cuts at the 14 TeV LHC. Negligible backgrounds are indicated by -. 



a(fb) 


Cuts Eq. Q + Q 


+ Eq. (0) 


■ pass n 

+ rij =2 


+ Eq. Q 


+subjet cuts 


Wh 


84 


27 


22 


3.3 


1.1 


Wjj 


1.4 x 10 6 


5.3 x 10 4 


4.4 x 10 4 


450 


1.2 


WW 


1.7 x 10 3 


41 


31 


0.91 


1.3 x 10~ 2 


wz 


470 


50 


42 


13 


1.2 x 10~ 2 


tq 


7.1 x 10 3 


170 


140 


0.28 


- 


tw 


3.3 x 10 3 


150 


76 


3.5 


3.3 x 10^ 2 


tbW 


2.3 x 10 5 


7.5 x 10 3 


1.5 x 10 3 


78 


- 


tt 


4.3 x 10 4 


1.4 x 10 3 


96 


26 


- 


S/B 


5.1 x 10~ 5 


4.4 x 10~ 4 


4.9 x 10" 4 


6.1 x 10~ 3 


0.85 


S/^S + B (30 fb" 1 ) 


0.36 


0.60 


0.57 


0.78 


3.8 


S/^S + B (100 fb" 1 ) 


0.66 


1.1 


1.0 


1.4 


7.0 



by noting that higher order QCD processes can generate a jet mass. For jets originating from 
massless partons, the average jet mass is approximately 



K) 



C^p 2 T R 2 , 



TT 



(3.3) 



where R is the jet radius and C depends on the initial partons. Hence, jets with lower px have 
on average lower masses than higher px jets. This effect can be seen in Fig. 0(b), where the 
second hardest jet is more (less) likely than the hardest jet to have a jet mass below (above) 
rrirj = 12 GeV. Based on these observations we apply the invariant mass cuts 

90 GeV < rrijj < 110 GeV and 8 GeV < rrij < 16 GeV. (3.4) 

The effect of the invariant mass cuts on the signal is shown in the second column of Table [j]. 

The irreducible backgrounds are the QCD background Wjj and the electroweak back- 
grounds WW/WZ — > Ivlj. Other contributing reducible backgrounds are tq — > iubj, 
tW -» £ub2j, tt -»■ £ubb2j, and tbW -»■ £u2b2j. The effects of the cuts in Eqs. Q, Q, 
and ( |3.4|) on the background cross sections are shown in the second and third columns of 
Table |]. Note that the tt processes are included in the tbW background and are shown only 
for reference. As can be seen, the invariant mass cuts greatly reduce all backgrounds. To 
further reduce the irreducible backgrounds, we require that the number of jets, n? ass , passing 
the cuts in Eq. Q3.1| ) to be equal to two. The effects of the nj ass = 2 cut on signal and 
background are shown in the fourth column of Table [l|. 

After the previous cuts, Wjj is still the dominant background. Equation ( |3.3| ) implies 
that the effect of the jet mass cut in Eq. (|3.4j) is to cause the QCD jets to strongly peak at 
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Figure 1: Invariant mass distributions of the signal Wh — > £v2j at the 14 TeV LHC for (a) the 
two hardest jets and (b) the individual jet invariant mass for the hardest (solid) and second hardest 
(dashed) jets. Cuts in Eq. (3.1) and (3.2) have been applied. 



low pt with a short tail into the the high px region. Whereas, since the signal jets originate 
from massive particles, the jet mass cut is expected have much less effect on the signal jet pt 
distributions. Hence, the signal jets are expected to have longer tails into the high px region 
than the background jets. In Fig. [2] we show the transverse momentum distribution of the 
(a) hardest and (b) second hardest jets that pass the cuts in Eqs. (jO|), (3.2), and (3.4) for 
both our Wh signal (solid) and Wjj background (dashed). As expected, the background jets 
peak at lower px and have shorter tails into the high px region than the respective signal jets. 
Hence, we apply the further px cuts on signal and background jets: 



Py(hard) > 100 GeV and ^(second hardest) > 50 GeV 



(3.5) 



The effect of these cuts on signal and background are shown in the fifth column of Table [l]. 

After all the above cuts, the signal is still difficult to observe over the background. With 
a luminosity of 30 fb _1 the significance at the 14 TeV LHC is still only 0.78er, improving only 
slightly to 1.4(7 at 100 fb _1 . As mentioned previously, our signal jets originate from color sin- 
glet massive particles while the background jets originate from colored massless partons. We 
now exploit those differences and perform a subjet analysis to increase the signal significance. 



3.1 Jet Substructure 



For the subjet analysis we follow the standard procedure given in Ref. |12]. First, for the two 
jets passing all previous cuts, the final step in the jet reconstruction is reversed to determine 
the two leading subjets: j\ and ji with m,j x > rrij 2 . Since the parent jet originates from a 
massive particle and the subjets from massless partons, there is expected to be a significant 
mass drop between the subjets and parent jet. Hence the subjet masses are required to meet 
the criteria rrij 1 < firrij, where rrij is the parent jet's mass and \x is a free parameter. Also, 
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Figure 2: Transverse momentum distribution of the two jets passing the cuts in Eqs. (3.1), (3.2), and 
(3.4) for (a) hardest and (b) second hardest jets at the 14 TeV LHC for both the Wh signal (solid) 
and Wjj background (dashed). 



for the decay rj — >■ cc the splitting should not to be too asymmetric 



AK ji,32 > Vent, 



in 



(3.6) 



j 



where AR 2 • = (eft 31 — (ft 32 ) 2 + (y 31 — y^ 2 ) 2 measures the separation of the subjets in the 
y — (ft plane and (ft 3 * is the azimuthal angle of subjet jj. If the subjets do not satisfy the mass 
drop and asymmetric splitting criteria, the parent jet is replaced by j\ and the procedure is 
repeated. As noted in Ref. [JDJ, for fi > l/\/3 the three-body decay r\ — > ccg will pass the 
mass drop criterion if the decay is in the Mercedes configuration in the rj rest frame. Following 
Refs. gO], Q we take fi = 0.67 and y cut = 0.12. 

Each jet in our signal is expected to result from the rj — >• 2c decay, where the two 
quarks will be collimated into a single jet. Hence the KT distance of the two subjets d = 

. 9 -9 

min(py ,py )AB% j 2 /R 2 ~ 0{m 2 ) pi|| . Figure y shows the log yd distributions for the 
(a) hardest and (b) second hardest jets with d measured in GeV 2 and after the cuts in 
Eqs. ( |3.1|) , (|3.2| ), and (|3.4| ) are applied. The Wh signal distribution is shown with solid lines 
and the Wjj background with dashed. Since we apply the jet mass cut in Eq. (|3.4| ), both 
the signal and background distributions are peaked near the same value of log yd. However, 
since the signal jets have a natural mass scale of m^ = 12 GeV, the signal distributions are 
peaked narrowly at log yd ~ log 12 ~ 2.5. Also, since the Wjj jets originate from massless 
partons, the background distribution has a significant tail in the low log yd region and is 
peaked at a slightly lower value than 2.5. We therefore apply the subjet cuts: 



2 < log Vd < 3. 



(3.7) 



Since our signal jets originate from color singlet particles, the two subjets within each jet 
are expected to form a color dipole. Hence, for our signal most of the radiation is expected 
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Figure 3: log \fd distributions at the 14 TeV LHC for (a) the hardest jet and (b) second hardest jet 
for both Wh signal (solid) and Wjj background (dashed). Cuts in Eq. (3.1), (3.2), and (3.4) have 
been applied and d is measured in units of GeV 2 . 
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Figure 4: Dipolarity distributions of the signal Wh — > tv2j at the 14 TeV LHC for (a) the hardest 
jet and (b) second hardest jet for both Wh signal (solid) and background (dashed). Cuts in Eqs. (3.1), 
(3.2), and (3.4) have been applied. 



to lie between the two subjets; whereas, the radiation in the background jets is expected to 
be more uniformly distributed since the jets originate from colored partons. The observable 
"dipolarity" [^] has been proposed to measure whether or not subjets form a color dipole: 

1 



V 



AR l v . 



E^» 



(3.8) 



where the sum runs over all calorimeter cells within a jet, p^ is the pt in a given calorimeter 
cell, and pTj is the jet's total pr- The distance between a given calorimeter cell and the line 
connecting the centers of the two subjets in the y — (f> plane is given by Ri. Fig. shows the 



dipolarity of the (a) hardest and (b) second hardest jets after cuts in Eq. (|3.1[) , (|3.2|) , and 
( |3.4[) . Again, the Wh signal is shown with solid lines and the Wjj background with dashed 



lines. For this distribution the dipolarity has been calculated at the hadronic level. As can 
be clearly seen, the distribution for the signal jets is peaked at low dipolarity and the Wjj 
distributions are more uniformly spread. We therefore apply the dipolarity cut 

V < 0.1 (3.9) 

The effects of the cuts of Eqs. (^^) and (|3.9[) are shown in the last column of Table |]. 
Including the subjet cuts improves the statistical significance of the signal to 3.8a at 30 fb _1 
and 7.0a at 100 fb _1 of integrated luminosity at the 14 TeV LHC. A statistical significance 
of 5a can be obtained with approximately 50 fb _1 of data. 

4. Conclusions 

The LHC is quickly excluding the remaining parameter ranges for the SM Higgs boson. 
Together with LEP direct detection limits, the only remaining available masses for a low mass 
SM Higgs are 117.5 - 118.5 GeV and 122.5 - 129 GeV |, |, ||. As these limits tighten, we 
are motivated to investigate other mechanisms of EWSB. Particularly interesting models are 
those that escape the LEP bounds and can have a Higgs boson with a mass around 100 GeV. 
This can be accomplished if the Higgs boson undergoes a cascade decay h—*2r}—> 4X, where 
j] is a pseudoscalar and X are light SM particles [|| . 

We have studied the observability of a Higgs boson in the "charming Higgs" model []12|| . In 
this model the Higgs decays via h — > 2-q — > Ac and can escape the LEP bounds. In particular, 
we considered the associated production of a Higgs boson with leptonically decaying W for 
masses rrih = 100 GeV and m v = 12 GeV. Since rrih 3> m^, the decay products of r\ — ^ cc are 
highly collimated and form single jets. Our signal then consists of a charged lepton, missing 
energy, and two jets originating from the two 77s from the Higgs decay. Such a signal is 
typically overwhelmed by the Wjj QCD background. By applying cuts to the invariant mass 
and transverse momentum of the jets, as well as requiring that exactly two jets satisfy our 
basic acceptance cuts, we can partially lift this signal above the standard model background 
at the LHC, but still only reach S/^S + B of 0.78cr at 30 ftr 1 and 1.4a at 100 fbr 1 at the 
14 TeV LHC. 

Our analysis can be greatly improved by applying jet substructure techniques. We require 
a significant mass drop from jet to subjet and fairly symmetric splitting of momentum between 
the two subjets. Additionally, we place new cuts on the KT distance between our subjets 
(which should be peaked around mi for our signal) and the jet dipolarity [|^] (because our 
signal comes from color singlets whereas our background is dominated by colored particles). 
Including these subjet cuts allows us to achieve a statistical significance of almost 4a at 30 
fb _1 and 7a at 100 fb _1 at the 14 TeV LHC. A 5a discovery can be achieved with 50 fb _1 of 
data. 

The procedures presented here are quite general and applicable to situations in which 
the Higgs boson undergoes the decay h — > 2rj — > 4j with m v <C rrih- Also, the branching 
ratio of h — > 4c may differ from one and the Higgs mass may be lager than 100 GeV. Simply 



scaling the signal rate with branching ratio, we estimate that with 100 fb _1 of data at the 14 
TeV LHC, branching ratios of BR(/i — > 277 — > 4c) = 0.35 and 0.64 are observable at the 3<r 
and 5a levels, respectively. Also, as the Higgs mass increases the dijet invariant mass cuts 
will need to be tuned and the signal rate will decrease slightly due to the higher final state 
invariant mass. However, the favored signal dijet mass is further from the Z-pole and the 
dijet invariant mass cuts will be more efficient in reducing the WW and WZ backgrounds. 
With the tuned dijet invariant mass cuts, the QCD background is expected to decrease at 
least as quickly as the signal rate. Hence, a slightly higher Higgs mass is not expected to 
significantly alter our conclusions. 
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